Bound State of Ultracold Fermions in the Spin Imbalanced Normal State 
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We study a polarized Fermi gas and demonstrate that Fano-Feshbach (FF) resonances lead to 
the pairing of fermions and holes into long living massive bosonic modes (bifermions and biholes), 
which can be viewed as signatures of a first order quantum phase transition. These modes, which 
are neither conventional Cooper pairs nor closed channel bosonic dimers, have a peculiar dispersion 
relation and appear already in the weak coupling limit, i.e. for large FF detunings. Our result is in 
accord with claims of Schunck et al. [cond-mat/0702066] that by stabilizing the normal state the 
population imbalance still allows effect of s-wave fermionic pairing. 
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Quantum phase transitions in polarized Fermi gases 
have been a subject of investigation for many years [1-5]. 
In Rcf . [1] , it was demonstrated that spin polarization of 
electrons in a superconductor leads to a first order quan- 
tum phase transition. An applied magnetic field induces 
a Zeeman splitting Ez in the energies of spin up and spin 
down electrons and imbalances the equilibrium densities 
of the two species. Eventually, for large enough Ez the 
superconducting ground state becomes a polarized nor- 
mal metal rather than a superconductor. 

Let us define the imbalance in the system as J = 
jy^jv" ' "^T an< i N± being the number of spin up and 
spin down fermions. To form an s-wave Cooper pair, 
two fermions have to be in different spin states. There- 
fore, for J = 1 the system has normal ground state. On 
the other hand, the ground state of a Fermi gas with 
no population imbalance (J = 0) is superfluid. Since 
no crossover between two states of different symmetry 
is possible, the system will remain normal even at zero 
temperature (T = 0) as long as J exceeds certain critical 
value J c . 

In superconductors, this translates into the existence 
of a critical magnetic field: in particular, the phase dia- 
gram in the presence of Zeeman (acting only on the spins) 
magnetic field [6] includes the BCS phase for E z < A, 
a cohesistence region A < Ez < \/2A, and the normal 
state for Ez > V^A, where A is the BCS superconduct- 
ing gap at zero temperature in the absence of Zeeman 
splitting. The elusive Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state [2, 3] appears in a narrow region below 
Ez = 1.508A: in this state, the mismatch of the spin 
up and spin down Fermi surfaces induces the electrons 
to form pairs with finite momentum. As a result, the 
superconducting gap becomes a periodic function of the 
space variables. 

Experimental studies of superfluidity have been per- 
formed in ultracold Fermi gases [7, 8] in the BEC-BCS 
crossover, tuning the scattering length through Fano- 
Feshbach (FF) resonances [9, 10]. Effect of the polariza- 
tion was achieved by creating imbalanced mixtures of ul- 
tracold 6 Li atoms in two different hyperfine states (which 



we will refer to as spin states til)- In the presence of the 
imbalance, anomalous density profiles [11] and the dis- 
appearance of quantum vortices [12] have been reported. 
While in metals orbital effects associated with the mag- 
netic field can obscure the effects of Zeeman splitting, in 
neutral atomic gases this problem does not exist. More- 
over, while in metals the relaxation due to spin orbit 
interaction tends to equilibrate the spin populations, in 
ultracold atomic gases the spin is a good quantum num- 
ber, and the initial imbalance in the system is conserved. 
In these systems, the difference between the numbers 
and Ni of majority and minority atoms can be encoded 
by the difference 2h in the chemical potential of the two 
spin states, which plays the same role as the Zeeman 
splitting for electrons in metals. 

In fact, the transition from superfluid to normal state 
in polarized Fermi gases (unlike the transition in an un- 
polarized gas) is an example of first order quantum phase 
transitions. Being discontinuous, first order quantum 
phase transitions appear to have no critical manifesta- 
tions. However, as it was pointed out in [13], the signa- 
tures of such phase transitions are the anomalies (massive 
modes) in the excitation spectrum. These modes lead to 
sharp features in the observables and can determine the 
collective quantum dynamics. Moreover, the study of 
such modes may provide valuable information about the 
relaxation processes. 

The identification of the massive modes and evaluation 
of their spectrum is a non-trivial model dependent task. 
The purpose of this Letter is to carry out this program 
for the normal state of resonantly coupled fermions with 
population imbalance. We prove that arbitrarily narrow 
FF resonances lead to the appearance of massive modes 
if J exceeds the critical value J c so that the normal state 
is stable. We identify two kinds of bosonic quasiparticles: 
bifermions and biholes, having spin zero, associated with 
the pairing of fermions and holes. These modes have 
different features both from the closed channel dimers of 
the BEC-BCS crossover and from conventional Cooper 
pairs. They have a strongly J-dependent gap (shown in 
Fig. 1), and a peculiar dispersion relation, their energy 
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FIG. 1: Excitation gaps fif,/(g = 0) and £lbh(q = 0) of 
bifermions and biholes in units of Fermi energy as functions 
of detuning 6 for fixed imbalance J = = 0.8 (panel 

A) and as functions of the imbalance J for 8 = 2.5cf (panel 

B) . The inset shows the behavior of the gaps for small im- 
balance, in the same units, together with the boundaries of 
the continuum (dashed lines). Here, the resonance width is 
7 = 0.1. The symmetry in the spectrum of bifermions and 
biholes (Qbf ~ ~^hf) is recovered for large detuning in panel 
A and for small imbalances in panel B. 



being a decreasing function of momenta. 

The FF resonant Fermi gas can be modeled by the 
following grandcanonical hamiltonian: 

H= fdr PA (-toV) I + tfe» (-ihV) b 

L a 

+ gR.pJdrJ2[^ 1 JfXfk], (1) 

<7\<T2 

4 = - V + 4*rh; e b = f- -2fi + 6 , 
2m 4m 

where a y is a Pauli matrix, H.p. denotes the Hermi- 
tian part of an operator, fl(r), [f a (r)] create [annihilate] 
fermionic atoms with mass m and spin a = ±1/2, while 
b'(r), [b(r)] create [annihilate] bosonic molecules; [x is the 
chemical potential. The bare detuning So depends on the 
high momentum cutoff pa — > oo. It is related to the po- 
sition of the FF resonance 6, which is independent of po, 
and is connected with the s-wave scattering length a s by: 
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The key dimensionless parameter is 7 
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ratio of the width of the FF resonance to the Fermi en- 
ergy ep (cf = \i at J — 0). The model (1) suggests that 
in the weak coupling limit (7 <C 1) that the system would 
be in the normal phase as long as 5 > S C (J) ~ 2ep (Fig. 
2). However, current experiments performed with wide 
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FIG. 2: (Color online) Schematic picture of the phase dia- 
gram at small coupling (for a complete discussion, see [5]). 
Arrows indicate the trajectories corresponding to the plots 
in Fig. 1. The FFLO phase lies in the red region at the 
superfluid/normal state boundary. 



resonances (7 > 1) have shown that the system remains 
in the normal state down to much smaller values of the 
detuning: at resonance (6 = 0), J c was found as small as 
0.7 [12]. Wc believe that our conclusions remain qualita- 
tively valid for wide resonances and at resonance, as long 
as the normal state is stable. 

In the case of narrow resonances (7 < 1), we can em- 
ploy the one loop approximation for the bosonic Green 
function T>(x±, X2), which amounts to the Dyson equation 
in Fig. 3, V° and G° being the non-interacting Green 
functions for the bosons and fermions correspondingly. 

Since T> is nothing but the two fermion scattering am- 
plitude, bound states and resonances of the system are 
found by identifying its poles. Solving Eq. c) of Fig. 3 
by means of Eqs. a) and b), after integration over fre- 
quencies, we find the equation for the dispersion relation 
f2(g) of excitations in the particle-particle channel at zero 
temperature 
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Here, h and the chemical potential /i are the Lagrange 
multipliers corresponding to the conserved quantities 
JVf — N± and JVf + iVj_. In our discussion, h and /i can 
be calculated by setting 7^0, since any correction to 
the zero coupling limit would translate into corrections 
to Eq. (3) of higher order in the small parameter 7 <C 1. 
Hence, given the imbalance J and the density of atoms 
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n, one can determine h and fi from 
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Equation (3) has two solutions fl b f(q) and Qbh(q) 
which correspond to the two kinds of quasiparticlcs. 
These excitations have gaps, which can be found from 
the equation for f2 
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The dependence of the gaps fib f and fij,/, on J and 5 
is plotted in Fig. 1. 

Clearly, these quasiparticles have spin zero, and obey 
bosonic statistics. However, they are quite different from 
conventional Cooper pairs in that their energy is a de- 
creasing function of momenta and their excitation gap 
strongly depends on the imbalance J of the system. 
Bifermions and biholes are also quite different from the 
closed channel bosonic dimers of the BEC-BCS crossover: 
in the limit 8 3> 2ep, such dimers are too energetic to play 
any role in the low energy physics of the Fermi gas. In 
contrast, bifermions and biholes can be observed experi- 
mentally even at large 8. Their excitation gap Qq ~ 2h 
can be by orders of magnitude smaller than the binding 
energy 8 of the bare bosonic dimers. 

In the limit J< 1, the quasiparticle spectrum can be 
evaluated analytically. Equation (3) becomes 
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with h = \J^f 



Here, A is the T = supcrfluid order 
parameter, which satisfies the gap equation 
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For 8 > 2e F , A ~ 8e F exp[-2 

The dispersion relation of these quasiparticles is non- 
trivial. For J <C 1, bifermions and biholes turn out to 
have opposite spectra. In this limit, as well as in the 
limit 8 3> e F (see Fig. 1), bifermions and biholes can be 
viewed as particles and antiparticles. For q <C J^/me Fl 
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FIG. 3: Diagrammatic representation of the Dyson equation 
for the particle-particle propagator, in the zero temperature 
formalism. 




the dispersion relation is given by: 



n bf (q) = -n bh ( q ) = n bf (o) 
n v (« = o) = 



l + 2(» b/ (0)/A) 2 2 2 

m b/ (o) Vpq 




(7) 



where vp = y/2e F /m. According to Eq. (7), which is our 
main analytic result, the excitation energy is a decreasing 
function of momenta, and the excitation gap depends on 
the coupling between the fermions through the superfluid 
order parameter A. 

In the one loop approximation, bifermions and biholes 
are undamped up to the critical value of momentum q c 
[14] , where their spectrum intersects a cut of the bosonic 
Green function (Fig. 4). For q > q c , the collective modes 
develop a nonzero imaginary part. 

The vanishing of Q b f(q = 0) and Q. b h(q = 0) at 
J = f — manifests the instability of the normal state 
with respect to formation of the homogeneous A for 
J < | A. . This corresponds to the first order phase tran- 
sition mentioned above. As to the instability toward in- 
homogeneous A(q) <~ e* qr , it appears at J c = 1-13^ [3], 
corresponding to a second order quantum phase transi- 
tion to the FFLO phase with broken translational sym- 
metry. 

The critical behavior near the second order transition 
between the normal and FFLO state at J < J c , as well 
as the collective excitations, definitely deserve theoreti- 
cal analysis. Here we can note just that the bifermion 
and bihole modes have little to do with these problems. 
Namely, the analysis of the solution of Eq. (3) on the 
physical and nonphysical sheets suggests that the FFLO 
transition manifests itself in the appearance of new dif- 
fusive branch rather than in softening of the mode (7). 
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FIG. 4: (Color online) The main plot shows the dispersion 
relations of bifermions and biholes, for J — J c (frequencies 
are in units of 2h, momenta in units of 2/i/Vf)- The value q c 
for which the spectrum intersects the cuts f2 = ±(2/i — VFq) 
(red lines in the main plot) is plotted in the inset as a function 
of Je^/A, in units of A/vf- 



In the context of disordered superconductors, the ex- 
istence of massive modes was discovered through mea- 
surements of the tunneling density of states [15], and the 
theory of the effect was constructed in Ref. [13]. We 
believe that bifermions and biholes in ultracold atomic 
gases should manifest themselves in measurements of 
radio-frequency spectroscopy which have detected super- 
fluidity at J = [8]. Possible additional experimental 
applications of those modes will be discussed elsewhere 
[16]. 

In conclusion, we have shown that FF resonances lead 
to collective modes in a Fermi gas with spin imbalance. 
As long as the imbalance is large enough to prevent the 
formation of inhomogeneous configurations, the physics 
is dominated by the appearance of bosonic quasiparticles 
(bifermions and biholes). These quasiparticles have finite 
gap, anomalous dispersion relation, and exist as long as 
the system is in the normal state. 

Added Note. After this paper was completed, experi- 
mental data were published by Shunck et al. [17], sug- 
gesting that s-wave pairing of atoms indeed takes place in 
the spin imbalanced normal phase of ultracold fermionic 
gases, as an explanation of anomalies in the population of 
imbalanced 6 Li atoms after applied radiofrequency. We 
find it likely that bifermion modes, which we discussed 
here provide a coherent interpretation of those anomalies. 
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